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ABSTRACT 

The main challenge for understanding the fuelling of supermassive black holes in ac¬ 
tive galactic nuclei is not to account for the source of fuel, but rather to explain its 
delivery from the boundaries of the black hole sphere of influence (10-100 pc) down to 
sub-parsec scales. In this work, we report on a series of numerical experiments aimed 
at exploring in further depth our model of “overlapping inflow events” as catalysts for 
rapid accretion, seeding a turbulent fleld in the infalling gas. We initially set a gaseous 
shell in non-equilibrium rotation around a supermassive black hole. After infall, the 
shell stalls in a disc-like structure. A second shell is then set in either co-rotation or 
counter-rotation with respect to the first and is let to impinge on the previously-formed 
disc. We And that combined turbulence and overlap significantly enhance accretion in 
counter-rotating inflows, while turbulence dominates for co-rotating inflows. The left¬ 
overs of overlapping inflows are warped nuclear discs, whose morphology depends on 
the relative orientation and angular momentum of the disc and the shell. Overlapping 
inflows leave observational signatures in the gas rotation curves. 

Key words: quasars: super massive black holes - galaxies: evolution - galaxies: active 
- black hole physics - hydrodynamics - turbulence. 


1 INTRODUCTION 

Understanding accretion onto super-massive black holes 
(SMBHs), from galactic scales of several kpc down to the 
horizon scale of few /i-pc, is critical for understanding their 
growth and determining their relation with the host galaxy 


(|Soltan 

1982| |Salucci et al. 

19991 |King 

2003| Shankar et al. 

2004IIM 

arconi et al.|2004||]V 

erloni|2004 

IHopkins et al.|2006| 


Raimundo Fabian|20Q9 ). Gravitational instabilities, either 
internal or triggered by mergers, are believed to be the main 
catalysts of SMBH growth over cosmic time, giving rise to 
the so called SMBH-galaxy co-evolut ion scenario ([Ferrarese 


Merritt |2000||Gebhardt et al.|2000||Ferrarese fc Ford|2Q0^ 

Kormendy Ho||2013| ). 

Once gas makes its way to the sub-parsec region, 
magneto-hydrodynamical processes, as discussed by |Balbus| 
& Hawley (1998), drive gas down to the innermost stable 
orbit of the SMBH. However, AGN accretion discs must 
be continuously or episodically replenished on these scales. 


* E-mail:jcarmona@sissa.it 


because at larger distances the discs are expected to frag¬ 
ment due to self-gravity (e.g. Shlosman et al.|[l990 Good- 
man||20b3l |Lodato||2012 ). For thin accretion discs, with as¬ 
pect ratio H/R ^ 10~^, gravitational instabilities start be¬ 
ing important at a distance of the order of Rqg — lO^Rs, 
where Rs is the Schwarzschild radius. For a black hole of 
Mbh = 10® M©, Rsg is equivalent to 0.01 pc. However, it 
is still a matter of debate to find suitable mechanisms that 
can make the gas cross the angular momentum barrier in 
the vicinity of the gravitational influence radius of the BH, 
Rinfl = GMbh/ct^ ~ 11 (Mbh/ 1 O®M 0 ) pc, where a ~ 100 
km s“^ is the velocity dispersion of the stellar component. 

Several mechanisms able to bridge this gap have been 


proposed in the literature (see Alexander & Hickox 


2012 ), 


which basically fall into two groups (Bellovary et al. 2013): 
either i) the gas falls almost radially in cold filaments from 
large distances leaving a sizeable amount of gas close to the 
SMBH (Bournaud et al. unu |Di Matteo et al.||2012| [Barai 


et al.||2012 Gaspari et al.||201^ Feng et al.||2014 ) or ii) the 

gas loses angular momentum by a series of processes such as 
gravitational instabilities, stellar feedback, turbulence, dis¬ 
ordered flows or the formation of nested bars ( Wada|[2004 


© 0000 RAS 


















































2 Carmona-Loaiza et al. 


odato fc Natarajan |2006l [Kawakatu fc Wada| 


Jogee| 


20Q8| [Hopkins fc Quataert|2010| ^oi et al.|2Q13| ). In partic- 

ular and interestingly for the purposes of this paper, |Hobbs| 
et ah (2011) have shown that supersonic turbulence is a 


key feature for enhancing the fuelling rate onto a SMBH 
from 100 pc scales down to the central parsec, because it 
creates convergent high density flows. Dense filaments are 
set into ballistic trajectories that travel through the ambi¬ 
ent gas, ultimately hitting the inner boundary. [Hobbs et al.| 
(2011) model their flows as initially unstable rotating spher¬ 


ical shells to capture the physics of this basic process. 

In a previous study (Carmona-Loaiza et al. |2014| Paper- 
I hereafter), we explored the hydrodynamics of a converg¬ 
ing non-turbulent spherical shell of gas impacting on a pre¬ 
existing stalled gaseous disc, adopting the initial condition 


proposed by Hobbs et al. (2011). This model, albeit highly 


idealised, was constructed to study how the interaction and 
collision of a converging flow onto a disc associated with a 
previous “failed” episode of accretion would enhance the ac¬ 
cretion rate onto a SMBH, to complement the [Hobbs et ah] 
(2011) model. The disc was not implanted in situ and con¬ 


structed from an equilibrium model, but generated from an 
infalling shell with a given angular momentum. The level of 
CO- or counter-rotation between the impacting fluid and the 
disc was, not surprisingly, found to play a key role in deter¬ 
mining the amount of gas dragged down to sub-pc scales: a 
major inflow of gas is triggered when the shell is counter¬ 
rotating relative to the disc. The interesting aspect is that 
not all of the gas is dragged inwards, but a remnant of nested 
rings, warped and inclined at various angles, is left. In con¬ 
trast with the usually proposed mechanisms, the overlap¬ 
ping inflow scenario provides a mechanism not of angular 
momentum redistribution through torques but rather of an¬ 
gular momentum cancellation via shocks. 

In this paper, we aim at studying the role of turbulence 
in misaligned overlapping inflows, and at highlighting poten¬ 
tial observational signatures of the fuelling mechanism pro¬ 
posed. We test the stability of the stalled nuclear disc against 
self-gravity in the Appendix as its fragmentation would 
reduce the hydrodynamical interaction between the flows. 
These tests allow us to select the maximum mass for which 
fragmentation is unimportant within our simple model. Our 
suite of hydrodynamical simulations is constructed consid¬ 
ering as main parameters; i) the initial angular momentum 
of the uniformly rotating shell, proportional to its rotation 
velocity r^rot; h) the relative orientation between the angular 
momenta of the disc and the shell, iii) the level of tur¬ 
bulence seeded in the fluid, measured by t’turb = {0,r’rot}; 
and iv) the self-gravity of the disc, related to the ratio 
m = Mdisc/TfeH. Actual astrophysical inflows are far more 
complex (see e.g. Wada 2004), however, to pursue such a 
task, one must be able to discern how gas behaves under 
changes of the most general parameters studied here. Tur¬ 
bulence plays the critical role of widening the angular mo¬ 
mentum distribution of the gas in both the disc and the shell, 
enlarging the cross section of interaction between them, and 
the outcome of this interaction cannot be predicted in ad¬ 
vance. 

The paper is organized as follows: In section we de¬ 
scribe the model and the computational method. In partic¬ 
ular we illustrate the implementation of the turbulent veloc¬ 
ity spectrum and of the time dependent artificial viscosity 


in the code. In section [S] we collect the results of all the sim¬ 
ulations to quantify the amount of mass crossing our inner 
boundary around the SMBH, and in section we describe 
the flow dynamics of the turbulent shells in the presence of 
stalled nuclear discs. We focus on the morphologies obtained 
and how they relate to the angular momentum content of 
the shell relative to the disc. A way of observing these relic 
structures through the analysis of the line of sight velocity is 
discussed in section and, finally, we draw the conclusions 
of the study in section Appendix]^ is devoted to testing 
the impact of self-gravity on the sub-pc inflow rate and on 
disc formation from single infalling shells. 


2 INITIAL CONDITIONS 


The environment in which our system evolves resembles 
that of the nuclear region of a galactic bulge. Specifically, a 
SMBH is located at the origin of an isothermal non-singular 
spherical stellar bulge described by a mass distribution 


M(r) = Mbh + 


Afcore {t j T core) ■) 
-Tfbulge ('^/^b ulge) 


T <C. T core 
'f' ^ 'f'corei 


( 1 ) 


with Mbh = lO^M©, Mcore = 2x lO^M©, Mbuige = lO^^M©, 
Tcore = 20 pc, Tbuige = 1 kpc. This mass distribution does 
not represent an active component in our simulations. It 
is just a static background potential into which the gas is 
immersed to explore the gas dynamics. While it has been 
shown that a stellar component can actually exert gravita¬ 
tional torques on the gas ( Hopkins &: Quataert|2010~ Chang 


et al.[[2007[ ), we here focus on a particular mechanism, not 
of angular momentum transport but instead of angular mo¬ 
mentum cancellation. The addition of a live potential (i.e. 
stellar particles) would introduce other angular momentum 
transfer processes which would make less clear the impact 
caused only by the interaction between the inflowing shell 
and the nuclear disc. 

We model the gas as a set of Asph — 10® SPH 
(Smoothed Particle Hydrodynamics) particles, in which half 
of them belong to the nuclear disc and the other half be¬ 
long to a shell surrounding it which is already unstable to 
inflow. For the simulations in which just a single shell is 
evolved, Asph — 5 x 10®. We introduce an accretion radius 
race = 1 pc surrounding the SMBH. Gas particles crossing 
this radius with low enough mechanical energy and angu¬ 
lar momentum are regarded as eventually accreted by the 
SMBH within a few orbital times and thus are removed from 
the simulation ( Bate et ^^[1995 ). 

The initial conditions for the overlapping inflows are 
very similar to those presented by [Carmona-Loaiza et al.[ 
(2014), this time giving the shell an additional turbulent ve¬ 
locity field superimposed on the rotational one in the way 


described by Hobbs et al. (2011). The initial conditions of 


each simulation consist of a disc-like structure which extends 
from the accretion radius race = 1 pc out to r ~ 20 pc, and 
a spherical shell with uniform density which extends from 
Hn = 30 pc to rout = 100 pc. Each shell is given a cylin¬ 
drical velocity profile = rrot = const., which can take 
the values rrot — 42,63 or 146kms“^ for the different cases 
treated here. In addition to the rotational velocity, the par¬ 
ticles forming the shell are given a turbulent kick spectrum. 
The mean turbulent velocity is not a free parameter and 
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is tied to the rotation velocity: I’turb = The turbulent 
velocity chosen in this way is such that it maximises the 
angular momentum distribution without overtaking the dy¬ 
namics (see Hobbs et al.|2011 ). 

The disc is always generated from a shell like the one 
just described having initially Vrot = '^^turb = 63kms“^. This 
single shell (i.e. without having a disc interior to it) evolves 
as it is not in equilibrium initially.Gas falls towards the cen¬ 
tre and angular momentum mixing due to shocks combined 
with the turbulent velocity field, leads to the formation of an 
extended disc (as opposed to the narrow ring in Carmona- 
Loaiza et ah 2014). After a time of t ~ 1.5Myr the disc 


reaches a stationary angular momentum distribution and 
the accretion rate becomes almost constant in time. At this 
time, we stop the simulation and surround the disc by a new 
uniform-density shell like the one from which the disc was 
born. The new shell, apart from having a different rotation 
and turbulent velocity, has its rotation axis tilted with re¬ 
spect to the axis of rotation of the disc. The latter defines 
the axis in all of our simulations, and so the disc initially 
lies in the x — y plane. We will refer to the disc in its initial 
condition as the primitive disc throughout this work, and its 
particles will be distinguished from shell particles. 

An isothermal equation of state has been assumed, 
P = CsP, with a constant temperature of T = lO^K. To focus 
on the hydrodynamical interaction between the two streams 
of gas under study, we regard the gas as non-self gravitating. 
Nevertheless we perform six additional simulations of single 
shells taking into account the self-gravity of the gas to asses 
its effects and the validity of the non-self-gravity approxi¬ 
mation. In the shell + disc simulations the mass of the shell 
is equal to the mass of the disc, Mgheii = Tfdisc- The units of 
mass, length and time that we will use for our computations 
are Mu = lO^^M©, Lu — I kpc, and Tu — 'o Myr respec¬ 
tively. The corresponding unit of velocity is Vu — 208kms“^ 
and the gravitational constant is automatically set G — 1. 
From now on we will continue to use these code units unless 
otherwise specified. 

The turbulent velocity spectrum was given as a Gaus¬ 


sian random field with a Kolmogorov power spectrum (Kol- 
mogorov]|1941 ), 


Pv{k) - k 


- 11/3 


( 2 ) 


where k is the wave number, since this is still the simplest de¬ 


scription that we have of a turbulent fluid (see however Fed- 
errath (2013) and references therein). The key assumption 
here is that the velocity field is homogeneous and solenoidal, 
so that we can define v in terms of a vector potential A such 
that F = V X A. As such, A is also described by a Gaussian 
random field with an associated power spectrum. 


PA{k) - k 


- 17/3 


( 3 ) 


To avoid a divergence in the variance of |A| at small /c, we 
introduce a cut-off /cmin, so that 


PAik) = Cik^ + klA 


,- 17/6 


( 4 ) 


where C is a normalization constant given to ensure that, 
given a box of volume 8 rgheii containing the shell, 


87rrfh,,ii 

(2^)3 


f 

Jkrr 


PA{k) k^ dk = -llturb. 


( 5 ) 


Physically the cut-off /cmin can be interpreted as the scale 
Amax = 27r//cmin ~ 'Tsheii, the largest scale on which the 
turbulence is likely to be driven and /cmax = 10^/cmin was 
introduced to cut the integral when convergence has been 
reached. 

To generate the statistical realisation of the turbu¬ 
lent velocity field we first sample the vector potential A 
in Fourier space using a periodic cubic grid of dimension 
256^. We draw the amplitudes of the components of Ak at 
each point {kx,ky,kz) from a Rayleigh distribution with a 
variance given by and assign phase angles that are 

uniformly distributed between 0 and 27r. We then Fourier- 
transform the vector potential to real space and take its curl 
to obtain the components of the velocity field. Finally, for 
estimating the turbulent velocity field at the position of each 
SPH particle, we use tricubic interpolation. The parameters 
of each simulation are given in Table 

It is worth saying that the precise form of the power 
spectrum is irrelevant for the scenario studied, as the im¬ 
portance of turbulence lies in widening the final distribu¬ 
tion of the gas. Nevertheless, the power spectra of supersonic 
solenoidal turbulence and subsonic turbulence have the same 
scaling relation with one important difference: while in sub¬ 
sonic turbulence the energy spectrum of the velocity field 
P[v] oc in supersonic (solenoidal) turbulence, it is not 

the velocity field that has the same scaling, but th e density- 
weighted velocity field, such that P[p^^^v] oc k~^^^ (Lighthill 


1955| |Henriksen||1991| |Fleck||1996| [Kritsuk et al.||2607| |Fed- 

erratli|2013| ). Since in our initial conditions the gaseous shell 
is homogeneous and the velocity field is solenoidal, the en¬ 
ergy spectrum to use is indeed indistinguishable from that 
predicted by Kolmogorov for incompressible homogeneous 
turbulence. 

We used the Lagrangian SPH code P- Gadget ( Springel] 
2005 Dolag et al.|2005 ) with a time dependent artificial vis¬ 
cosity (AV) as described in Morris & Monaghan ( 1997| ) (see 
also [Valdarnini 2011). In its standard formulation, the ar¬ 
tificial viscosity a^isc always acts when modelling turbulent 
flows, as a good deal of convergent streams exist. However, 
when applied to turbulent (supersonic) fluids, the expecta¬ 
tion is that the standard viscosity ends up damping the tur¬ 
bulent velocity spectrum unphysically fast. Thus, we include 
a time dependent artificial viscosity aviso (t ): a viscosity pa¬ 
rameter that rises to non-negligible values only where and 
when it is needed (i.e. at a shock), decreasing to almost zero 
anywhere else. In this scheme, each particle is assigned its 
own AV parameter, ai which is evolved according to the 
differential equation: 


dai 

dt 


Oii OLu 


+ Si. 


(6) 


Thus, ai always decays to amin with e-folding time r^, while 
the source term Si will make its value rise rapidly when 
approaching a shock. For our isothermal equation of state 
P = CsP, the decay time-scale is given by (Price & Federrath 


2010 ): 
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Table 1. Parameters of the simulations: 4 simulations of single 
shells with the same specific angular momentum content (given by 
Vrot) and varying mass (Mgheii) made in order to asses the impact 
of self-gravity; 6 simulations of single shells of the same Mgheii 
with varying Vrot and turbulent velocity (nturb)^ fo quantify the 
relevance of turbulence for dragging mass towards the central 
parsec; 12 more simulations of overlapping inflow events in near 
CO- and counter-rotation = 60° and 150° respectively) with 

and without a turbulent velocity field to isolate the boost to the 
inflow rates produced by the overlapping of the inflows. 


run 

r’rot 

Vtnvh 

Otiit 

Mgheii/M bh 

V03-M2 

0.3 

0.3 

- 

0.01 

V03-5M2 

0.3 

0.3 

- 

0.05 

V03-M1 

0.3 

0.3 

- 

0.1 

V03-M0 

0.3 

0.3 

- 

1.0 

v02-nt 

0.2 

- 

- 

0.1 

v03-nt 

0.3 

- 

- 

0.1 

v07-nt 

0.7 

- 

- 

0.1 

v02 

0.2 

0.2 

- 

0.1 

v03 

0.3 

0.3 

- 

0.1 

v07 

0.7 

0.7 

- 

0.1 

v02t060-nt 

0.2 

- 

60° 

0.025 

v02tl50-nt 

0.2 

- 

150° 

0.025 

v03t060-nt 

0.3 

- 

60° 

0.025 

v03tl50-nt 

0.3 

- 

150° 

0.025 

v07t060-nt 

0.7 

- 

60° 

0.025 

v07tl50-nt 

0.7 

- 

150° 

0.025 

v02t060 

0.2 

0.2 

60° 

0.025 

v02tl50 

0.2 

0.2 

150° 

0.025 

v03t060 

0.3 

0.3 

60° 

0.025 

v03tl50 

0.3 

0.3 

150° 

0.025 

v07t060 

0.7 

0.7 

60° 

0.025 

v07tl50 

0.7 

0.7 

150° 

0.025 


Ti = 


hi 

(JCi ’ 


( 7 ) 


where cr is a non-dimensional free parameter which deter¬ 
mines over how many information-crossing times the viscos¬ 
ity decays. For all of our simulations we have taken cr = 0.1. 
The source term is implemented as 


Si = So fi max(0, -V • Vi), (8) 

where So is another free parameter that we set to So = 
1.0. A maximum value for the viscosity parameter is also 
introduced to keep the viscosity bounded between amin = 
0.02 and CTmax = 1.0. 

The (hydrodynamical) smoothing, /i, and (gravita¬ 
tional) softening, £, are both adaptive and coupled in such 
a way that e — h. We impose a minimum softening de¬ 
pendent on the mass of the shell simulated (in all cases 
£min <0.1 pc). In Appendix[^we select the maximum mass 
for which fragmentation of the disc is unimportant within 
our simple model. For the rest of the simulations in which 
calculations were done neglecting self-gravity the gravita¬ 
tional softening plays no role and the smoothing length is 
fully adaptive with no threshold /imin being imposed. 



^rot/^disc,0 


Figure 1. Fraction of the mass (in units of Mg^eii) in the cen¬ 
tral parsec as a function of the initial rotational velocity for all of 
the simulations in this work and those of Paper-I. Colours blue 
and green refer to shells falling in co- and counter-rotation with 
the disc; circles and squares refer to the absence or presence of 
a turbulent kick in the initial conditions; empty symbols refer to 
shells falling in the absence of a nuclear disc. Thus, from bottom 
to top, the data points joined by lines represent the following: i) 
single shells with no turbulence (empty circles; < 10“^; shown 
also in the inset), ii) co-rotating overlapping inflows with no tur¬ 
bulence (green circles; < 10“^; shown also in the inset), iii) co¬ 
rotating overlapping turbulent inflows (green squares), iv) tur¬ 
bulent single shells (empty squares), v) counter-rotating over¬ 
lapping inflows with no turbulence (blue cirlces), vi) counter¬ 
rotating overlapping turbulent inflows (blue squares). It is evi¬ 
dent that the presence of a nuclear disc signiflcantly enhances the 
mass flowing into the central parsec when the ~ 10 — 100 pc scale 
inflow is counter-rotating with it. See text for details. 


Table ^ lists the three sets of simulations, the first for 
a single converging shell, the second describing the interac¬ 
tion of the shell with a pre-existing disc with no turbulence, 
and the third with turbulence included. The results of the 
first set with self-gravity are illustrated in the Appendix and 
motivate the choice of using as reference value for the mass 
of the shell and the disc, Mgheii = Mdisc = 0.025Mbh. 


3 TURBULENT VERSUS IMPACT INFLOW 
MODE 

Figure ^summarises our main result using as key indicator 
the “accreted” mass Mace- We plot the fraction of Mgheii 
that crossed the central parsec after t = 0.3, versus the ini¬ 
tial rotation velocity of the shell, Urot (expressed in units of 
the rotational velocity of the shell which formed the stalled 
disc, Udisc,o = 0.3) for the entire data set (See table [^. From 
the figure we can infer two different modes of angular mo¬ 
mentum cancellation and redistribution: One driven mainly 
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by turbulence -referred to as turbulent mode, and one driven 
by the impact of the shell with the disc -referred to as im¬ 
pact mode. 

The data set comprises cases in which: 

• (I) a single non-turbulent shell falls inwards onto the 
SMBH (black empty circles); 

• (II) a single turbulent shell falls inwards onto the SMBH 
as in [Hobbs et al. ( 2011| (black, empty squares); 

• (HLa) a non-turbulent shell interacts with the stalled 
disc in near co-rotation (green, hlled circles); 

• (HLb) a turbulent shell interacts with the stalled disc in 
near co-rotation (green, filled squares); 

• (IV. a) a non-turbulent shell interacts with the stalled 
disc in near counter-rotation (blue, filled circles); 

• (IV. b) a turbulent shell interacts with the stalled disc in 
near counter-rotation (blu, hlled squares). 

The sets of parameters that deliver the smallest fraction 
of gas towards the central parsec are those with no turbu¬ 
lence and no interaction (I). Right above these -although not 
noticeable in the Figure, non turbulent shells in co-rotation 
with a stalled disc (HI.a) deliver just < 10“^Msheii of gas 
to the central parsec. The yield is increased to as much as 
15% of Msheii for turbulent single shells (H). The interaction 
with the stalled disc in near co-rotation lowers the accreted 
mass of the turbulent shells (cases HLb). This was not ex¬ 
pected and can be explained noting that the denser, higher 
angular momentum fluid elements present in the disc drag 
and accelerate the more rarefied fluid elements. These shell 
elements with lower angular momentum are those that hap¬ 
pen to interact soonest with the disc, and this prevents them 
from reaching the accretion radius. 

The two most efficient sets of parameters for driving gas 
to the central parsec are those in which the inflow is driven 
by the impact-mode (cases IV.a and IV.b). In the cases 
(IV.a), up to 1/3 of the shell’s gas is poured into the central 
parsec, and the yield is raised up to Maccc ~ l/2Msheii in 
the cases (IV.b). Even for the shells with the highest angular 
momentum (urot/'^’disc,o = 2.33), the impact-mode of accre¬ 
tion drags a non-negligible Mace ~ O.lMgheii amount of gas 
into the central parsec. Even if their high angular momen¬ 
tum would intrinsically make these shells deliver a negligible 
amount of mass into the central parsec, the interaction with 
the disc raises this number to levels comparable with those 
produced by lower angular momentum shells. 

These considerations hold true if the masses involved 
are low enough for the nuclear discs not to disrupt (see Ap¬ 
pendix]^, as more massive clumpy discs would be expected 
to evolve into stellar discs and consequently would not ex¬ 
ert the required hydrodynamical torques for triggering sub- 
pc inflows from the 100 pc scales. It is the presence of the 
gaseous disc that causes shell fluid elements impacting on it 
to cancel their angular momentum, ultimately driving gas 
towards the central parsec. Turbulence magnifies this effect 


by i) creating converging flows, as described by Hobbs et al. 


(2011), and by ii) widening the section of interaction be¬ 
tween the shell and the disc. 


4 RELIC STRUCTURES 

In this section, we describe the properties of the relic struc¬ 
tures that form after the violent shell-disc interaction. 


The key role of turbulence is that of broadening the 
disc formed by the first infalling shell. This is because tur¬ 
bulence widens the initial angular momentum distribution 
(Hobbs et al. 2011), causing gas to circularise in lower as well 
as in higher angular momentum orbits relative to the mean 
circularization radius. Some of this gas crosses the accretion 
radius and some remains at the periphery forming a wide 
disc. When a second shell interacts with the relic disc, the 
collision proceeds roughly as described in Paper-I: the shell 
impacts the disc and the angular momentum of both disc 
and shell change under the condition of angular momentum 
conservation of the shell+disc system. This time however, 
the interaction takes place over a wider radial range. Thus, 
while the structures studied in Paper-I lead to the forma¬ 
tion of tilted nested rings, in this work the presence of an 
initial turbulent kick leads to the formation of continuous 
warped discs surrounded by a more extended and dynami¬ 
cally decoupled disc. As we have seen, the inflow rates are 
also enhanced. 


4.1 Counter-rotating flows 

Eigure shows the complex structure that forms after the 
hydrodynamical interaction, in the counter-rotating inter¬ 
action case^ As the gas of the shell rains onto the disc, it 
shock-mixes with the fluid particles of the disc. The almost 
head-on encounter leads to an important cancellation of an¬ 
gular momentum, dragging a significant amount of mass 
within the accretion radius. Eluid particles that still pre¬ 
serve enough angular momentum form an extended, nested 
warped disc comprising both shell and disc particles. In none 
of these cases do particles inside 5 pc rotate at the initial 
tilt angle (^tiit = 0). The right panels of Eigurej^show the 
angle ^ as a function of spherical radius r for comparing the 
strength of the warps. 

• In the run v02tl50, the inflowing shell has a smaller an¬ 
gular momentum than that of the disc (Jsheii = 2/3 Jdisc)- 
Two decoupled warped discs are formed. One inside the 5 
central parsecs, and the other extending from 5 to 15 parsecs 
approximately. Both have lost memory of the initial angular 
momentum distribution. 

• In run v03tl50, the inflowing shell has the same angular 
momentum as the primitive disc (Jsheii = Jdisc)- At first 
sight (Eigure|^ left central panel), it seems that there are 
three disconnected discs. However, in fact only two decou¬ 
pled discs are present (right central panel of Eigure]^: one 
warped disc extending from 10 pc inwards and an outer 
non-warped lighter disc extending from 15 to 30 pc. One 
would naively expect to have all the gas of the shell mixed 
with all the gas of the disc generating a single massive disc, 
shrunk and with its angular momentum making an angle 
0-1/2 X 150° with the axis. However, since the interac¬ 
tion starts taking place as soon as the first gas from the shell 
impacts on the disc, the formation of the inner warped disc 
occurs before the rest of the shell has time to make its way 
down to the stalled disc. This late fraction of the shell still 
interacts with the rest of the disc, changing its orientation 
and shrinking its orbits but without producing significant 


^ The snapshots of our simulated discs were done using SPLASH 
( Price|2Q07| ) 
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inflow. At this time, because of the shrinkage of the primi¬ 
tive disc, particles of the shell that were on the high tail of 
the angular momentum distribution will not be able to inter¬ 
act any more, thus remaining in orbit as an outermost less 
dense disc with the same orientation the shell had initially 
((9tiit = 150°). 

• In run v07tl50, the inflowing shell has a higher angular 
momentum than that of the stalled disc (Jsheii = 7/3 Jdisc)- 
Two nested decoupled discs are formed with only the inner 
one being warped. The inner disc extends from the accre¬ 
tion radius out to 15 pc, while the outer non-warped disc 
extends from 15pc outwards (Lower panels of Figure [^. The 
large kinetic mismatch between the two flows leaves behind 
a decoupled disc, with a non trivial distribution of angular 
momentum which might reflect also how we set our initial 
conditions. 


4 . 1.1 Co-rotating inflows 

When the infalling shell is nearly in near co-rotation with 
the stalled disc, there is no significant cancellation of angu¬ 
lar momentum. The interaction leads to the formation of a 
single warped disc-like structure in which 0 changes orien¬ 
tation smoothly with radius, i.e. a larger scale warped disc 
forms as illustrated in Figure For the cases in which the 
inflow has lower and higher angular momentum than the 
stalled disc (i.e. runs v02t060 and v07t060), the warp looks 
particularly clean, whereas for the cases in which the inflow 
and the primitive disc have the same angular momentum 
(run v03t060) the disc is more disturbed and the warp more 
structured. Note that 0 is not changing monotonically with 
distance r for any of the cases plotted. 

For run v03t060, a large number of particles is concen¬ 
trated just below and just above the circularization radius 
(^circ = 16 pc) with rotation axes spanning a wide range of 
angles at the same radius (see bottom middle panel of Fig¬ 
ure]^. These are particles whose rotation axis is either tilted 
by ^ ~ 10°, almost parallel with the z axis of the primitive 
disc, or tilted by ^ ~ 50°, almost aligned with the mean 
angular momentum of the inflow. This suggests that the in¬ 
teraction has not yet finished, even if most of the accretion 
has already taken place. Gas at these radii will continue to 
shock-mix until a state of equilibrium is reached in which 
the particles at each radius rotate around the same axis. 

Runs v02t060 and v07t060 are smoother. At each radius 
the particles rotate around a single axis (see Figure left 
and right lower panels). From the inside to the outside the 
rotation axes in both cases span almost AO ~ 60°. 

As a general trend, gas at the circularization radius of 
the primitive, stalled disc has its angular momentum aligned 
with the z axis while gas at the circularization radius of the 
shell has its angular momentum tilted with respect to the z 
axis by keeping memory of its original orientation. 


5 OBSERVATIONAL SIGNATURES 

A comparison between the rotation profiles of a single disc 
and of a warped disc formed after the interaction of two 
nested inflows is carried out in Figures and where we 
show a random sample of 10 % of the particles for two se¬ 
lected simulations: a single disc generated by an infalling 


turbulent gaseous shell (run v03) and the end state of a shell 
that fell onto a stalled disc with lower angular momentum 
and in near counter-rotation (run vrt02tl50). The left and 
right panels in each figure correspond to the gas observed 
from two perpendicular lines of sight (LOS) that coincide 
with the primitive disc being edge-on and face-on respec¬ 
tively. In the top panels the LOS velocity of each particle 
is plotted against its radial distance. In addition, a red line 
showing the velocity profile, vk oc cos(z), expected for 

a Keplerian disc, is superimposed, where the angle z = 0 
for an edge-on disc and i = 90° for a face-on disc. From 
these, we immediately note that the single turbulent disc 
(Figure]^ has a rotation profile consistent with an almost 
flat Keplerian disc, whereas the composite case of a shell in 
counter-rotation impacting on the disc (Figure shows a 
very different profile, inconsistent with Keplerian rotation. 
In the lower panels, LOS velocity maps are shown with red 
and blue colours indicating receding and approaching gas. 
In the single disc case (Figure]^, a very ordered rotation is 
observed: gas approaching from the left and receding to the 
right when viewing the disc edge-on, and gas with a very 
tiny almost constant velocity when the disc is face-on. In¬ 
stead, when observing the remnant of the overlap (Figure 
1^, the ordered pattern just described is not present and it 
is clear that two kinematically decoupled bodies of gas are 
present. 

The comparisons just made must be taken only 
as an indication of the kind of signatures expected. 
To obtain more realistic signatures, one should take 
into account a live stellar potential which would cer¬ 
tainly exert torques on the gas, thus contributing to 
the final shape of the warp. An additional interest¬ 
ing problem is that of exploring the possible forma¬ 
tion of stellar warped discs, which must be tackled 
including in-situ star formation. 


6 DISCUSSION 


Previous works have examined the collision of a gas cloud 


with a pre-existing nuclear gas disc around a black hole (Alig 
et al.|2Q13 ), the collision of two gas clouds orbiting a central 
black hole ( Hobbs Nayakshin||2009 ) and the impact of a 
prolate gas cloud impinging and engulfing a central black 
hole ( Lucas et al.|20I3' ) by means of idealised hydrodynam- 
ical simulations with the particular aim of explaining the 
origin of the disc of young stars observed at the Galactic 
Gentre around Sgr. A* (e.g. Genzel et al.||20I(5 ). The com¬ 
mon feature of all of these experiments is the collision of gas 
flows with different angular momenta and the subsequent 
formation of warped and/or nested discs. If the Galactic 
Gentre represents a standard, these observations comple¬ 
mented by theoretical modelling suggest that galactic nuclei 
are sites where collisions between gas clouds and residual 
discs take place. The model which we have presented here 
and in Paper-I can be regarded as related to these processes. 
It has a higher degree of symmetry than those but also gen¬ 
erates warped and nested rings. It seems plausible then, that 
a relation between the mass inflow and the angular momen¬ 
tum difference between two gaseous components (disc-cloud, 
cloud-cloud) can be obtained, even if in an approximate way. 
This could be relevant for more accurate sub-grid models 
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Figure 2. Warped discs formed when an inflow of gas, here modelled as a rotating gaseous turbulent shell, impacts on a stalled disc 
in near counter-rotation = 150°). The left column shows the generic outcome of these events. The central column shows a 

zoom into the innermost warped disc that forms (note that the disc view has been changed here to highlight better the degree of disc 
warping). The right column shows 10% of the particles selected randomly. For each particle we record its distance r (in units of pc) 
from the SMBH and its tilt angle 0 between the speciflc angular momentum of the particle and the 2 : axis (representing the rotation axis 
of the primitive disc). Red particles come from the shell while blue particles come from the primitive disc. We see that the inner 
and outer disc are always in counter-rotation. From top to bottom, the rows correspond to the cases v02tl50, v03tl50 and v07tl50. The 
Figure illustrates the high degree of warping of the structure. 


of AGN feedback in cosmological simulations (e.g. Gaspari 
et al.|2013| |. 


Strong circumstantial evidence indicates that warps are 
common in AGN accretion discs. For instance, maser emis¬ 


sion in the nuclear regions of NGG 4258 (Herrnstein et al. 
2005), Gircinus ( Greenhill et al.|2003 ), and four of the seven 


galaxies examined by Kuo et al. (2011) has revealed the pres¬ 


ence of warped gaseous discs at 0.1-1 pc scales. Addition¬ 
ally, warped thin discs have been invoked by |Nayakshin fc] 
Guadra (2005) as the obscuring mechanism behind the uni¬ 


fication paradigm of AGN. They argue that warped clumpy 
discs provide the easiest way to obscure the central engine 
of AGN. 


All of this suggests that a connection between warped 
discs and AGN fuelling is not just accidental. However, 
we need to keep in mind that, in order to estab¬ 
lish a connection between the warping of the disc 
and AGN fuelling, a realistic treatment of the stellar 
potential must be included, as gravitational torques 
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Figure 3. Warped discs formed when an inflow of gas, here modelled as a rotating gaseous turbulent shell, impacts a primitive, stalled 
disc in near co-rotation = 60°). The upper row shows the generic outcome of these events: a large-scale warped disc which 

changes its plane of rotation from the inside to the outside. The lower row shows 10% of the particles selected randomly. For each 
particle we record its distance r (in units of pc) from the SMBH and its tilt angle 6 between the speciflc angular momentum of the 
particle and the 2 ; axis (representing the rotation axis of the primitive disc). Red particles come from the shell while blue particles 
come from the primitive, stalled disc. Note that these structures are always continuous discs changing their plane of rotation by as much 
as ~ 60°. From left to right, the columns correspond to the cases v02t060, v03t060 and v07t060. Note that for the case v07t060 it would 
seem that the inner part of the disc is detached from the outer part, however the r — 6 plot shows otherwise. In all cases there is a 
continuous mix of shell and primitive disc particles. 


exerted by it could potentially modify the structure 
of the warps and the AGN fuelling rate. 


7 CONCLUSIONS 

The presence of a stalled gas disc around an SMBH on scales 
between ~ 1 — 30 pc, where the angular momentum bar¬ 
rier is found, can enhance the inflow of matter that sub¬ 
sequently reaches these regions. By interacting via shocks 
with the stalled disc, angular momentum cancellation takes 
place. In this way, the stalled disc acts as an agent or “ob¬ 
stacle” that modifies the free flow of any gas reaching it. 
Misalignment between the angular momenta of the arriving 
gas and that of the disc is preferred for an efficient cancella¬ 
tion of angular momentum. Such a situation is represented 
in our suite of simulations by the shell reaching the disc in 
counter-rotation, which greatly enhances inflow (see Figure 
[^. An even greater inflow is obtained in cases in which the 
shell is initially given a turbulent velocity field. 

The interactions described in this work generate nested 
and kinematically detached warped discs and establish a 
natural connection between the angular momentum cancel¬ 
lation mechanism and the warping of the discs. Co-rot at ing 


inflows give origin to smoother warps that are nonetheless 
interesting to study. Such warped discs on ~ 10 pc scales, 
show line-of-sight velocity profiles that depart significantly 
from Keplerian rotation, as illustrated in Figure These 
warped structures, on parsec scales, arise from the purely 
hydrodynamical interaction between a misaligned inflow and 
a stalled gas disc. Neither SMBH feeding nor disc warping 
in this scenario invoke the action of gravitational or viscous 
torques (see e.g. Tremaine &: Davis|2Q14 ). More realistic sim¬ 
ulations of e.g. cloud-cloud collision, motivated by state of 
the art simulation of galaxies under formation and evolution 
will asses whether the “cancellation of angular momentum 
by a stalled gaseous disc” is an important catalyst of SMBH 
feeding and AGN fuelling. At the same time, more real¬ 
istic simulations including a live stellar potential will 
quantify the relative impact on the disc’s shape of 
the overlapping inflows scenario with respect to the 
gravitational torques generated by a non-spherical 
stellar bulge. 
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Figure 4. Single disc viewed edge-on (left) and face-on (right). 
In the top panels the line of sight (LOS) velocity of each particle 
is plotted in black against its distance from the SMBH. Keplerian 
LOS velocity profiles are superimposed in red. Note that in both 
cases the velocity profiles of the particles are consistent with Kep¬ 
lerian rotation: r’LOS oc cos(z), where i = 0 when the disc is 

edge-on and z = 90 when it is face-on. In the lower panels, we see 
the LOS velocity map. Bluer colours represent gas approaching 
the observer while redder colours represent gas which is receding. 
In the edge-on case, the LOS velocity goes from being positive 
on the left side to negative on the right side (with an absolute 
maximum value of ~ 500 kms“^), while in the face-on case the 
LOS velocity is very small. 



Figure 5. Warped discs that form after the shell-disc interaction 
for the case vrt02tiltI50 (the shell which fell with lower angu¬ 
lar momentum than the disc and in near counter-rotation). The 
lines of sight are the same as in Figure]^ In the top panels the 
line of sight (LOS) velocity of each particle is plotted in black 
against its distance from the SMBH. Keplerian LOS velocity pro¬ 
files are superimposed in red. Note that in both cases the velocity 
profiles of the particles are not consistent with a flat disc in Ke¬ 
plerian rotation: z;los oc cos(z), where i = 0 when the disc 

is edge-on and i = 90° when it is face-on. In the lower panels, we 
see the LOS velocity map. Bluer colours represent gas approach¬ 
ing the observer while redder colours represent gas receding from 
her/him. In both of the bottom panels it is clear that in the cen¬ 
tral region the gas is kinematically detached from the rest. 
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APPENDIX A: SELF-GRAVITY 


In this Appendix we explore the role of self-gravity by vary¬ 
ing the mass of the gaseous inflowing shell which will then 
settle into the “stalled” disc. For very low shell masses as 
compared to the mass of the SMBH, m = Msheii/AfeH << 1, 
we expect that the inflow will proceed cleanly as if the fluid 
were non self-gravitating. With increasing m, we expect the 
formation of clumps due to local gravitational instabilities. 
By doing this exercise, we determine the mass threshold be¬ 
low which self-gravity is negligible, and at the same time we 
assess the role of self-gravity in altering the values of the 
inflow rate onto the SMBH in our simulations. 

In a massive gaseous disc, the parameter that tilts the 
balance towards fragmentation or stability is the cooling 
time. Indeed, for a disc to be stable against self-gravity, the 


Toomre (1964) parameter. 


Q 


Cs K 

ttGE’ 


(Al) 


where k and S are the epicyclic frequency and the sur¬ 
face density respectively, must be greater than unity, with 
marginal stability occurring at Q = 1. When the cooling 
time-scale of the gas, tcooi, is greater than the dynamical 
time-scale of the disc, tdyn, an initially hot disc (Q > 1) will 
cool slowly until the marginal instability stage is reached 
(Q = 1). During this process, the disc has had enough time 
to respond, with local gravitational perturbations forming 
into a global spiral structure that in turn heats the gas 
through shocks and pdV work, preventing the disc from go¬ 
ing into the fragmenting regime (Q < 1). Thus, the interplay 
between local self-gravity (governed by tdyn) and the cooling 
of the gas (governed by tcooi) will dictate the evolution of the 
disc. In the scenario just described (tcooi > tdyn), the spiral 
structure acts as a thermostat, always maintaining the gas 
at marginal stability. On the other hand, if tdyn > tcooi, as 
the disc cools faster than the time which it needs to respond 
and develop a spiral structure, the disc breaks into clumps 
(for a more detailed discussion of these processes see e.g. 


|Lodato|2007l|2012 |. 


Gold constant-temperature discs, having an idealized 
tcooi = 0, are thus prone to fragmentation. In this way, the 
discs in our simulations, with a temperature T = 10^A, 
stand on the worst-case scenario in which there is no way 
of halting fragmentation once the disc becomes unstable. 
In a constant-temperature thin disc with a Keplerian rota¬ 
tional profile (k ~ ^(^)), Q depends on the mass of the 
disc only through its surface density, F. If we were to in¬ 
clude a more realistic cooling treatment, more massive discs 
would still be able to survive fragmentation, making it eas¬ 
ier to feed the black hole through overlapping inflows and 
thus strengthening our conclusions. The Toomre parameter 
calculated for a non self-gravitating disc is shown in Figure 
|A1[ The minimum of the Toomre parameter, Qmin, occurs at 
the circularization radius (rdrc — 16 pc) due to the greater 
concentration of mass there. 

We performed 4 simulations of single turbulent inflows 
in which the total mass in the gaseous shell was varied, tak¬ 
ing into consideration the self-gravity of the shell. The ini¬ 
tial conditions were the same as those described in section 
in these specific cases with an initial rotation velocity 
of Vrot = '^^turb = 0.3 and a mass Mgheii = s x Mbh, for 
s = 0.01, 0.05, 0.1, 1.0. 


Al The development of clumps following the 
infall of a single shell 

Glumps developed in all of the simulations of this section 
with the exception of the lowest mass run m = 0.01. In 
the left-panel of Figure |A2| we show clump formation for 
the case m = 0.1 at the time of disc formation (t = 0.3). 
The clumps develop within a narrow ring at a circulariza¬ 
tion radius, rdrc — 16 pc (as could be predicted by Figure 
Al). The clumps become gravitationally bound and revolve 
around the SMBH in quasi-Keplerian orbits. As they form at 
the circularization radius due to the concentration of mass 
at that radius, their rotation around the SMBH is quite 
stable, so that they barely interact with each other. When 
they do interact, they merge into larger and denser clumps, 
but no scattering occurs that can bring mass towards the 
SMBH. Tidal tails develop among them but these are not 
strong enough to disrupt the dynamics of the disc outside 
their vicinities. In the right-panel of Figure [A^ we show the 
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r [pc] 

Figure Al. The Toomre parameter calculated for a non self- 
gravitating thin disc in Keplerian rotation around the SMBH 
of mass Mdisc = O.IMbh- In this case, the surface density 
S(r) ~ 7 .- 372 ^ density p(r) ~ r~^ and the total scale height 
H(r) = 0.02 r. With this numbers at hand, the Toomre param¬ 
eter can be estimated as Q ~ p(r)] ~ A^MBn/^disc? 

with k ~ 0.1. The Toomre parameter is roughly constant within 
the inner 10 parsecs and reaches a minimum Q << 1 at the cir¬ 
cularization radius rdrc — 16 pc. Note that the same shape Q{r) 
shifts upwards or downwards depending on M^isc- 


smooth disc that forms when m = 0.01 corresponding to the 
case of negligible self-gravity. 

A2 Inflow Rates 

We now compare the specific inflow rates, M/Mgheii, to¬ 
ward the central parsec, produced by varying the mass of 
the inflow. By doing this we can have a grasp of how rele¬ 
vant self-gravity is in either helping or preventing gas from 
falling from scales of 10 — 100 pc down to the central par¬ 
sec. The specific inflow rates for the four cases are compared 
in Figure |A3| It appears that the specific accretion rates 
are independent of the mass of the inflow when the masses 
are low enough, Mgheii < 0.05Mbh- Even in the cases with 
Mgheii = 0.05 — 0.1 Mbh, in which self-gravity starts playing 
an important dynamical role by disrupting the inflow into 
clumps at Tcirc, the accretion rate during the most active 
phase, 0.03 < t < 0.2, remains unaffected by self-gravity. 
This means that we can safely model inflows of single shells 
as being non-self-gravitating whenever Mgheii ^ O.IMbh, as 
far as the sub-pc inflow rate is concerned. This threshold 
must change when dealing with overlapping inflows, as we 
must have a space-filling gaseous disc for the second inflow to 
interact with. Therefore, when treating overlapping inflows, 
we must have Mtot = Mdisc + Mgheii ^ 0.05Mbh. 
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Figure A2. Projected log(E/Msheii) in the xy plane when a disc form {t = 0.3 expressed in time units Tu = ^ Mry). The left-panel 
shows self-bound clumps in the Mghell = 0.1 Mbh model where self-gravity is important. The right-panel shows the disc in model 
^shell = 0-01 ^BH where self-gravity is unimportant. 



Figure A3. Particle accretion rates as a function of dimen¬ 
sionless time for the cases in which the mass of the inflowing 
shell is Mghell = Mbh (red dotted), O.IMbh (blue-dashed), 
Mghell = O.OSMbh (green dashed) and Mghell = O.OIMbh 
(magenta-dashed). For comparison, the particle accretion rate for 
a case in which the non self-gravity approximation was used is also 
shown (black solid line). The overlap of the lines makes it nearly 
impossible to distinguish between them at relevant levels of ac¬ 
cretion (t < 0.18; shaded region), which correctly illustrates that 
particle accretion rates are independent of the mass of the shell. 
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